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ABSTRACT 



Aims. In this paper we study the main spectral and temporal properties of gamma-ray bursts (GRBs) observed by Fermi/GBM. We 
investigate these key properties of GRBs in the rest-frame of the progenitor and test for possible intra-parameter correlations to better 
understand the intrinsic nature of these events. 

Methods. Our sample comprises 32 GRBs with measured redshift that were observed by GBM until August 2010. 28 of them belong 
to the long-duration population and 4 events were classified as short/hard bursts. For all of these events we derive, where possible, the 
intrinsic peak energy in the vF v spectrum (Ep fsX ), the duration in the rest-frame, defined as the time in which 90% of the burst fluence 
was observed (r^est) and the isotropic equivalent bolometric energy (E iso ). 

Results. The distribution of E pfest has mean and median values of 1 . 1 MeV and 750 ke V, respectively. A log-normal fit to the sample 
of long bursts peaks at ~ 800 keV. No high-£ p population is found but the distribution is biased against low E p values. We find the 
lowest possible E p that GBM can recover to be ss 15 keV. The Tgg^ distribution of long GRBs peaks at ~ 10 s. The distribution of 
E lso has mean and median values of 8.9 x 10 52 erg and 8.2 x 10 52 erg, respectively. We confirm the tight correlation between £p_ res t 
and E iso (Amati relation) and the one between £ p . rest and the 1-s peak luminosity (L p ) (Yonetoku relation). Additionally, we observe a 
parameter reconstruction effect, i.e. the low-energy power law index a gets softer when E p is located at the lower end of the detector 
energy range. Moreover, we do not find any significant cosmic evolution of neither E p iest nor r 9 o, re st- 

Key words. Gamma-ray burst: general 



1. Introduction 

Gamma-ray Bursts (GRB) are the most luminous flashes of y- 
rays known to humankind. It is generally believed that they orig- 
inate from a compact source with highly relativistic collimated 
outflows (T > 100). 

A large fraction of our knowledge of the prompt 
emission comes from the Burst and Transient Explorer 
(BATSE, Mee gan et all 1 19921) onboard the Compton Gamma- 
Ray Observatory (CGRO, 1991-2000). Unfortunately, only a 
handful of BATSE bursts had a measured redshift because the 
BATSE error boxes were too large and thus, follow-up observa- 
tions with X-ray and optical instrumentation was very limited. 
Moreover, the first afterglow was only detected in 1997, already 
near the end of the BATSE mission. The lack of distance mea- 
surements led to a focus of GRB studies in the observer frame 
without redshift corrections. Due to the cosmological origin of 
GRBs, such a correction is likely to be necessary to understand 
the intrinsic nature of these events. 

W ith the two dedicated satellite s, Beppo-SAX dBoella et alj 
1 19971) and Swift dGehrels et al.ll2004l) . the situation has changed 



and afterglow and host galaxy spectroscopy has provided red- 
shifts for more than 200 events by now. Unfortunately, the rela- 
tively narrow energy band of Beppo-5AX (0. 1 keV - 300 keV) 
and Swift/B AT (15 keV - 150 keV) limits the constraints on 
the prompt emission spectrum because the peak energy, £ p , 
in the vF v spectrum of GRBs, can only be determined for 
low E p values and i s often unconstrained dButler et al.1 [2007: 
ISakamoto et~aT1l2009l) . 

In this work we will take advantage of the broad energy cov- 
erage of the Fermi/GBM (8 keV - 40 MeV) to study the primary 
spectral and temporal properties, such as E p , Tgo, the time in- 
terval in which 90% of the burst fluence has been observed, and 
Ei so , the isotropic equivalent bolometric energy, in the rest-frame 
of the progenitors of 32 GRBs with measured redshift. 



2. GRB sample and analysis 

2.1. Instrumentation 

The Gamma-Ray Burst Monitor (GBM) is one of the in- 
struments onboard the Fermi Gamma-Ray Space Telescope 
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(lAtwood etal.ll2009t) launched on June 11, 2008. Specifically 
designed for GRB studies, GBM observes the whole unocculted 
sky with a total of 12 thallium-activated sodium iodide (Nal(Tl)) 
scintillation detectors covering the energy range from 8 keV 
to 1 MeV and two bismuth germanate scintillation detectors 
BGO) sensitive to energies between 150 keV and 40 MeV 
Mee ean et alJ l2009). Thus, GBM offers a unprecedented view 
of GRBs, covering more then 3 decades in energy. 

2.2. Burst selection 

The selection criterion for our sample is solely based on the 
redshift determination. We form a sample of 32 bursts detected 
by GBM up to October 16th, 2010, with known redshift (deter- 
mined either spectroscopically or photometrically). 

Our sample contains 4 short and 28 long GRBs. The red- 
shift distribution of the GBM GRBs is shown in Fig. Q] together 
with a histogram of all 239 GRBs with redshift determinations 
to dat^B A two-tailed Kolmogorov-Smirnov (KS) test between 
the full sample and the GBM-only sample shows that the two 
distributions are very similar (P = 84 %). In conclusion, the 
GBM-only sample is representative of the full GRB sample with 
redshift. 
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Fig. 1. Redshift distribution in % of GBM GRBs (blue solid line) 
compared to all 239 GRBs with measured redshift to date (red 
dashed line). Both samples contain long and short bursts. 



2.3. Analysis of the GBM Data 

For the E p determination, CSPEC data dMeegan et al.l2009l) with 
a time resolution of 1.024 s (4.096 s pre-trigger) were used. For 
the short GR Bs, i.e. those with Tgo < 2 s, time-tagged event 
(TTE, Meegan et alJl2009l) data were used with a fine time res- 
olution of 64 ms and the same channel boundaries as CSPEC 
data. 

Detectors with source angles greater than 60° and those 
occulted by the spacecraft or solar panels were discarded. A 
maximum of three Nal detectors were used for each E p de- 
termination. In four cases (GRB 090929A, GRB 090904B, 

1 www . mpe . mpg . de/ " j cg/grbgen . html 



GRB 090618, and GRB 081007) only one Nal detector could 
be used due to the reasons mentioned above. Where possible, 
both BGO detectors were included in the analysis if they were 
not occulted by the satellite during the prompt emission. Even 
if there was no apparent signal in the BGOs, they can help to 
determine an upper limit of the GRB signal in the spectra. 

The spectral analysis was performed with the software pack- 
age RMFIT0 (version 3.3rc8) and the GBM Response Matrices 
vl.8. To account for the changing orientation of the source with 
respect to the detectors caused by the slew of the spacecraft, the 
detector response matrices (DRM) were generated for every 2 
degrees on the sky. 

For each burst we fitted for every energy channel a low-order 
polynomial to a user defined background interval before and af- 
ter the prompt emission and interpolated this fit across the source 
interval. 

Three model fits were applied to all time bins with a signal- 
to-noise (S/N) of at least 3.5 above the background model: a sin- 
gle power-law (PL), a power law function with an e xponential 
high e nergy cutoff (COMP) and the Band function (Band et al. 
1993a). For three GRBs (GRB 090424, GRB 090618, and 
GRB 090926A) an effective area correction was applied to the 
BGO with respect to the Nal detectors to account for systemat- 
ics which dominate the statistical errors due to the brightness of 
these events. The best model fit is the fu nction which provides 
the best Castor C-stat value (Cash 1979). An improvement by 
AC-stat= 10 for every degree of freedom is required. The profile 
of the Cash statistics was used to estimate the l<x asymmetric 
error. 

The values obtained with this analysis method may be super- 
seded by the GBM spectral catalogue released by the GBM team 
(Goldstein et al., in preparation). 

3. The intrinsic properties 

3.1. The peak energy (E p ) 
3.1.1. Instrumental bias 

A first important issue which needs to be addressed is whether 
the bursts observed by GBM are drawn from the same distribu- 
tion as the bursts which were observed by BATSE. Because of 
the broader sensitivity of GBM to higher energies, there could be 
a significant deviation towards higher Ep values. To answer this 
question, we use the BATSE cataloguqf] from which we extract 
the E p values which were obtained from the time averaged spec- 
tra (fluence spectral fits). Ignoring the power-law (PL), Gaussian 
Log and Smoothly Broken Power Law (SBPL) fits, we elected to 
the COMP model or Band function for the purpose of this test. 
We require AE p /E p < 0.4 and the low-energy power-law index, 
a has an absolute error cr a < 0.4. GRBs which did not fulfill 
these criteria were rejected from the sample. The Band function 
was always preferred over the COMP model in cases in which 
the high-energy power-law index /3 was constrained (erg < 0.4). 
Both long and short GRBs were included in this analysis. The so 
obtained E p distribution shown in Fig. [2] 

The same selection cut was applied to the 2-year GBM spec- 
tral catalogue. The red histogram in Fig. [2] shows this distri- 
bution. Both distributions peak at * 170 keV and show the 
same standard deviation. A KS test reveals that the difference 

2 RMFIT for GBM and LAT analysis was devel- 
oped by the GBM Team and is publicly available at 
f ermi . gsf c . nasa . gov/ssc/data/analysis/. 

3 www.batse.msfc.nasa.gov/ goldstein/ 
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log E p [keV] 

Fig. 2. E p values of 1367 BATSE/CGRO GRBs with a log- 
normal fit (black dotted histogram and line), 375 GBM GRBs 
and log-normal fit (red dashed histogram and line) and 30 GBM 
GRBs with redshift measurement (blue solid histogram). A KS 
test suggests that all three samples were drawn from the same 
distribution. 



between the two samples is not statistically meaningful (P = 
18%). This means that the two histogr ams are drawn from the 
same distributi on, in agreeme nt with iNava et al .] (12010b . We 
note that, iBissaldi et al. I J2oTTh showed that the E p distribution 
of some GBM -GRBs exten ds to higher energies compared to 
BATSE (Kaneko et al. 2006). However, this is of no surprise as 
Amati et al. (2002) demonstrated that bursts with higher fluence, 
i.e. hi gher Ej so , have, on a verage, higher E p (see also Sect l4. 1ft . 
Since Bissaldi et a D (1201 lh only use bright GRBs, it is to be ex- 
pected that their E p distribution is shifted to higher energies. We 
conclude that GBM, although being sensitive up to 40 Me V, does 
not find a previ ously undiscove r ed pop ulation of high-Zip GRBs, 
consistent with lHarris & Shareldl998l) . 

The KS test was then applied to both the E p distributions be- 
tween the whole sample of GBM bursts and the 30 GBM bursts 
with measured redshift (for GRB 0905 19A and GRB 080928 a 
PL model fits the data best) and also to the E p distribution of 
BATSE bursts and the E p distribution of GBM bursts with mea- 
sured redshift. In neither case the differences were statistically 
meaningful (P = 24% and P = 20%, respectively). Thus, all 
3 histograms are very likely drawn from the same distribution. 
In conclusion, the sample of GBM GRBs with measured red- 
shift presented here is representative for the whole population of 
GRBs which were ever observed by BATSE and GBM. 

However, it should be stressed that GBM cannot measure 
E p values which are lower than a certain limiting threshold. It 
is w ell known that E n value s of GRBs can go as low as a few 
keV. iPelangeon et alj d2008l) for example find E p values as low 
as ~ 2 keV in GRBs that were obser ved by the High E nergy 
Transient Explorer 2 (HETE-2, see e.g. Barraud et allll2003l and 
references therein). These low energetic events have been classi- 
fied as X-ray flashes (XRF) or X-ray rich bursts (XRB) (see e.g. 
iHeise & in 'tZandl 12001: Sa kamoto et alJl20 05). However, it is 
very likely that XRFs and XRBs are nothing else than weak and 
long GRBs (see lKippen et aljr2004 and references therein). 



The borderline E p value is obviously located somewhere 
near the low-energy sensitivity of the Nals, which has yet to be 
determined. Thus, in order to determine a potential bias in the 
£p,rest distribution shown in Fig. [3] it is important to understand 
and quantify the limits of the GBM to measure E p (be it either 
from the COMP model or Band function). 

For this purpose, we created a set of simulated bursts with 
different initial spectral and temporal starting values. We input 
the source lifetime (ts, 1 s, 5 s, 10 s, 100 s) and the photon flux 
(/) in the 10 keV to 1 MeV range (1, 3 and 10 ph cirT 2 s -1 ). 
For the simulation the Band function was chosen as the photon 
model with varying E p (15,17, 25, 50, 100 keV) but fixed a = 
-0.8 and ft - -2.4. We simulate these bursts overlaid on real 
background data by using detector Nal 7 of GRB O9O926A0 
This results in 60 different burst models. For each model, we 
created 1000 bursts to account for Poissonian noise. This results 
in 60000 spectra, each of which was then fitted with the Band 
function using the detector response matrix (DRM) of detector 
Nal 7 created for the location of GRB 090926A. 

After the fitting procedure, we reject those bursts which have 
AE p /E p > 0.3 and cr a < 0.4. These rejected spectra are then 
defined as unconstrained. We did not apply this criterion to the 
high-energy power law index ft. A spectral fit that has a con- 
strained £p and a but an unconstrained ft is simply considered a 
COMP model. In TableQ]we report the mean and standard devia- 
tion of the output spectral parameters, E p and a, of the simulated 
bursts. 

The conclusions of this exercise are: 

1. GBM can recover E p values only as low as w 15 keV. The 
fact that the observed E p distribution is indeed b iased is in 
clear contradiction to e.g. iBrainerd et al ] (120001) . They ar- 
gued that the observed E p distribution (by BATSE in this 
case) is actually the intrinsic one. 

2. the input E p parameter can be recovered from the simulated 
spectra within the 2 <x errors for almost all simulated flux 
levels and source lifetimes. 

3. GRBs with low fluxes, low E p and short t s are more likely 
to be rejected than GRBs with higher fluxes, higher E p or 
longer t s 

4. the low-energy power law index a tends to get softer, i.e. 



to have lower values, for low E„ and low 



last point is particularly n o tewort hy. ICrider et al 



fluxes . Thi s 
(H997I). 



iLlovd-Ronning & Petrosianl (120021) and later Kaneko et al. 
(2006) found a significant correlation between E p and a in 
the time-resolved spectra of several GRBs. Supported by our 
simulations, we point out the possibility that a parameter 
reconstruction effect is at work in addition to any intrinsic 
correlation between Ey and a. This effect has already been 
brough t for ward bv iPreece et al.l (ll 998b : [Llovd & Petrosianl 
(2000) and ILlovd-Ronning & Petrosianl ( 120021) . In short, it 
depends on how quickly a can reach its asymptotic value 
and how close E p is located to the low-energy limit of the in- 
strument's energy bandpass. The closer E p is situated to the 
detector's sensitivity limit, the fewer is the number of pho- 
tons in the low-energy portion of the spectrum. This makes it 
increasingly difficult to determine the asymptotic value of a. 
Instead, a more negative value of the low-energy power law 
index will be measured which is what is observed here. 



4 The choice of Nal 7 and the choice for this specific GRB is com- 
pletely arbitrary. We could have chosen any other detector that observed 
any other real GRB. 
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Table 1. Mean and standard deviation of the output spectral pa- 
rameters of the simulated bursts. 
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This Table shows the input source lifetime, t s (in s), E p j n (in 



Notes, 

keV), photon flux, / (in ph cm 



s ), the recovered E„, £p, ou t ( m keV), 



the recovered low-energy PL index, a, and the percentage of bursts 
which were rejected due to the selection cut, rej (in %). The input 
low-energy power law index was -0.8. Cases for which all simulated 
bursts were rejected are not reported in this Table (for details see text). 



3.1.2. The intrinsic E p distribution 

In Fig. [3] we present a histogram of E pjest of our sample of 
GBM GRBs where a E p measurement was possible. The dis- 
tribution of all bursts has mean and median values of 1 . 1 MeV 
and 750 keV, respectively. A log-normal fit to the sample of long 
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Fig. 3. Rest-frame distribution of E p for 4 short (hatched his- 
togram) and 26 long (empty histogram) GBM-GRBs. For two 
GRBs (GRB 0905 19A, GRB 080928) a simple PL model fits the 
data best, thus not providing an E p value. The complete distribu- 
tion has mean and median values of 1.1 MeV and 750 keV, re- 
spectively. A log-normal fit to the sample of long GRBs (dashed 
line) peaks at ~ 800 keV and has a FWHM = 0.93 in log-space. 



bursts peaks at ~ 800 keV. In the canonical scenario of GRB jets, 
turbulent magnetic fields build up behind the internal shock, and 
electrons produce a synchrotron power law spectrum. The typi- 
cal rest-frame frequency v rest in the internal shock dissipation is 
v rest °c 0.2 Lgf r n 'MeV dZhang & Meszarosll2002l) . where L 52 
is the luminosity in 10 52 erg/s and the dissipation radius in 
units of 10 13 cm. Our measurement of the rest-frame peak en- 
ergy therefore leads to a constraint of L 1 ^ r~ 3 ' w 0.8 (in the 

above units). 

Recentlv. lCollazzi et all d2011l) reported that the width of the 
£p,rest distribution must be close to zero with the peak value lo- 
cated close to the rest-mass energy of electrons at 5 1 1 keV. This 
effectively implies that all GRBs must be thermostated by some 
unknown physical mechanism. We tested this claim by deter- 
mining the width of the E p KSl distribution (cteji+z) = 0.48) to- 
gether with all the individual errors that add to the uncertainty 
of the E p measurem ent in log-space (for details refer to eq. 1 in 
ICollazzi etalll201ll) : 



Total 



2 2 2 2 

: °"Poisson + °"Det + °"choice + °~Def ■ 



(1) 



cchoice describes the different choices made by different ana- 
lysts. Since all the bursts in this paper were analyzed consistently 
(selection of the time interval, energy range, etc.), we can set 
cchoice = 0. o"D e t is the error which results in not k nowing the de- 
tector response perfectly. However, according to ICollazzi et al.L 
cr Det is negligibly small, meaning that the calibrations of the de- 
tectors are usually well understood. <x Def describes the differ- 
ences that are obtained in E p when usi ng different photon mod- 
els. As it wa s show n already by e.g. iBand et aD (Il993bl) and 
iKaneko et alJ ( l2006h the C OMP model res ults in higher E p val- 
ues than the Band function. ICollazzi et alJ use cr Mo( i e i = 0.12 to 
account for this difference. Here, this value is adopted as <x De f . 
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The Poissonian errors are cr Po i sson 
0.16 



0.10 in log-space. Thus 



cr Total = V0.12 2 +0.1 2 

Inserting the just found values in eq. 8 in Collazz ret al.1 gives 



= cr 



E p (l+z) ' 



= 0.45. 



(2) 



However, a zero width in the E p distribution is synonymous 
with cr £pint = 0. The conclusion is that the £p jres t distribution 
does not have a zero width. This finding i s in conflict w ith the 
implications and conclusions discussed in Colla zzi et al.1 



3.2. GRB duration (Tgo) 

For determining the duration of a GRB, we applied the method 
first introduced by Kou vehotou et al. (1993). They defined the 
burst duration as the time in which 90% of the burst counts 
is collected (Tgo). Here, we adopted the same definition; how- 
ever, the burst's fluence was used instead of the counts. The 
T^o value depends highly on the detector a nd on the energy in - 
terval in which it is determined (see e.g. Biss aldi et al.l 120 lTh . 
Additionally, since we are interested in durations in the rest- 
frame of the GRB, it is not sufficient to simply account for the 
time dilation due to cosmic expansion by dividing the measured 
durations by (1 + z). The energy band in which the Tgo is de- 
termined needs to be redshift corrected as well. We determine 
the burst duration in fluence space in the rest-frame energy in- 
terval from 50 keV to 300 keV, i.e. in the observer frame energy 
interval from 50/( 1 + z) keV to 300/( 1 + z) ke V. 




12 3 

logT 9arest [s] 

Fig. 4. 79o jrest distribution of 4 short (red hatched histogram) and 
28 long (blue histogram) GBM GRBs in the redshift corrected 
energy interval from 50 keV to 300 keV. The log-normal fit to 
the distribution (dashed line) of the long bursts peaks at ~ 12 s 
and has a FWHM = 1.1 in log-space. 



3.2.1. The intrinsic Tg Q distribution 

In Fig|4]we present the rest-frame distribution of Tgo of our sam- 
ple of 32 GBM GRBs. The number of short bursts is still too 
small to unambiguously recover a bimodal distinction of short 
and long events in the rest-frame. Henceforth, a distinction is 
made between the short and long class of GRBs using an ad- 
hoc definition of 79o, res t < 2 s for short bursts and Tgo, rest > 2 s 
for long bursts, although we do have neither observational nor 
physical evidence to support such a distinction in the GRB sam- 
ple presented here. All GRBs in our sample that are defined 
as short in the observer-frame, also remain short with the here 
adopted definition. We note that GRB 1008 16A is peculiar in 
that it is classified as a short burst by GBM data, but it has a 
Tgo = 2.9+0.6 in the 15 keV - 350 keV energy r ange in Swift data 
with a low-level emission out to about T0+ 100 dMarkwardt et al.l 
120101) . 

The mean value of the whole Tgo jest distribution is 32.4 s and 
the median value is 10.8 s. 

3.3. Isotropic energy (E iso ) 
Eiso is calculated by 
And 2 L 

E'lSO — Z S y, (3) 

1 + z 

where d^ is the luminosity distance and S Y the fluence in 
the 1/(1 + z) keV to 10/(1 + z) MeV frame. We determine S y 
using the energy flux provided by the best-fit spectral parameters 
and multiplying it with the total time interval over which the 
fit was performed. Since we performed the fit for time intervals 
where the count rate exceeded a S/N ratio of 3.5, it happened 
that some time intervals of some bursts were not included in the 
fit (e.g. phases of quiescence where the count rate dropped back 



to the background level). These time intervals were not used to 
calculate the fluence. The 5 y distribution is shown in Fig|5] The 
median value of the fluence distribution is 1.6 x 10~ 5 erg cm -2 
and the mean value is 5.9 X 10~ 5 erg crrT 2 . A log-normal fit to 
the data peaks at 2.2 x 10~ 5 erg crrr 2 . 



7F 
6 \ 
5 'r 




log S [erg cm 2 ] 

Fig. 5. Energy fluence distribution of 32 GBM GRBs determined 
in the energy range between 1 /(l + z) keV to 10/(1 + z) MeV 
range. The complete distribution has mean and median values 
of 1.6 x 10~ 5 erg crrr 2 and 5.9 x 10~ 5 erg ctrT 2 . A log-normal fit 
(dashed line) to the long GRBs peaks at 2.2 x 10~ 5 ergcirT 2 and 
has a FWHM = 2.2 in log-space. 
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The Ziiso distribution is shown in Fig|6] The distribution for 
the long bursts has a median and mean value of 1.2x 10 53 erg and 
1 .4 x 10 53 erg, respectively. Short bursts, on the other hand, have 



significantly lower values of 2.9 x 10 M erg and 4.0 x 10 51 erg, 
respectively. A log-normal fit to the long bursts reveals a central 
value of 10 53 1 erg. Because our sample is dominated by long 
GRBs a log-normal fit to the whole distribution results in an es- 
sentially unchanged peak value (10 53 erg). 




50 51 52 53 54 55 56 
log E iso [erg] 



Fig. 6. Ei S0 distribution for long (blue hatched histogram) and 
short (red hatched histogram). For long bursts the distribution 
has mean and median values of 1.4 x 10 53 erg and 1.2 x 10 53 erg, 
respectively. A log-normal fit to the long bursts (dashed line) 
peaks at 1.2 x 10 53 erg with a FWHM = 1.9 in log-space. Short 
bursts have a mean and median value of 4 x 10 51 erg and 2.9 x 
10 51 erg, respectively. A log-normal fit to the combined sample 
(dash-dotted line) peaks at 9.8 x 10 52 erg with a FWHM = 2.2 
in log-space. 



4. Correlations 

4.1. Amati relation 



Amat i et aD d2002l) first showed that there is a tight correlation 
between £p, res t and E lso (the isotropic equivalent bolometric en- 
ergy determined in the energy range between 1 keV to 10 MeV). 
This relation is now known as the 'Amati relation". In FigJ7]we 
show the Amati relation for the 30 GBM GRBs with measured 
£p,rest and Ziiso. While there is an evident correlation between 
these two quantities (Spearman's rank correlation of p — 0.74 
with a chance probability of 1.7 x 10~ 5 ) the extrinsic scatter of 
the long GRBs is larg er by a factor of ~ 2 in log-space com- 
pared to lAmatil (1201 Oh Also, the best fit to our data is shifted to 
slightly larger £p ires t values. The best fit power-law index to the 
long GRBs of our sample is 0.5 2 + 0.06 which is in agreement 
with the indi ces obtained by e.g. lAmatH d20 1 Ot) : [Ghirlanda et al.l 
(2009| 1201 Oh. As has been shown by other authors in the past 
(see e.g. lAmati et al]f2008t Ghirlanda et al. l l2009HAmatill2010h 
short bursts do not follow the relation, being situated well out- 
side the 2 cr scatter around the best-fit. This is true also for the 
power-law fit derived here (see FigJT} except for GRB 1008 16A. 



However, as already stated above this burst may actually fall in 
an intermediate o r hybrid class of short GRBs with extend ed 
emission (see e.g. lNorris & Bonnell 20061: IZhang et al.ll2009l) . 




52 53 
log E iso [erg] 



55 



Fig. 7. Amati relation for 4 short (red open triangles) and 26 
long (blue open diamonds) GBM GRBs. Also shown is the best 
power-law fit to the data (black solid line) and the extrinsic 2<x 
sc atter (black d ashed lines) with the best power-law fit published 
by lAmatil d2010l) (red dash-dotted line) with the 2cr scatter (red 
long dashed line). 



4.2. Yonetoku relation 



lYonetokuet 10(12004 found a tight correlation between the rest 
frame peak energy in the vF v spectrum E pTest and the 1-s peak 
luminosity (L p ) in GRBs (so called Yonetoku relation). The peak 
luminosity is calculated with 



(4) 



where d\ is the luminosity distance and F p the 1 s peak en- 
ergy flux. We determine F p in the energy range between 30 keV 
and 10 MeV in the rest frame of the GRB (lYonetoku et al.l2 004). 
We use the GBM peak spectral catalogue (Goldstein et al. 201 1) 
to determine the best fit spectral parameters of the brightest 1 s 
bin. We then determined the energy flux in the rest-frame of the 
GRB by integrating 



f>10000/(l+z) 

F= No E d>(£) dE, 

J30/(l+z) 



(5) 



where iVo is the normalization of the spectrum (in photons 
cm 2 s _1 ), E the energy and 0(E) the form of the spectrum (either 
Band or COMP). In order to determine the error on the energy 
flux, we calculate 1000 flux values for each GRB by varying 
the input parameters, No, E p , a and ft according to the lcr error 
on each parameter. We then use the median and the mean abso- 
lute deviation (MAD) around the median of the 1000 simulated 
bursts to determine the peak energy flux and the error on the 
latter, respectively. 



6 



D. Gruber et al.: Rest-frame properties of 32 gamma-ray bursts observed by the Fermi Gamma-Ray Burst Monitor 



We present this relation for 30 GBM GRBs in Fig.[8]and find 
a best-fit power law index of 0.58 ± 0.08. The Spearman's rank 
correlation gives p - 0.7 with a chance pro bability of 2.3 x 10~ 5 . 
Our findings are in goo d ag reement with Yonetoku et al.l (|2004) ; 
iGhirlanda et al. (2009) and Ghirl anda et al.l (l2010t) . 




51 52 53 
log L p [erg/s] 



55 



Fig. 8. Yonetoku relation for 4 short (red open triangles) and 26 
long (blue open diamonds) GBM GRBs. Also shown is the best 
power-law fit to the data (black dashed line). 



4.3. rarest vs redshift 

The rest-fra me Tgg is plotted as a function of redshift in Fig. [9] 
Contrary to Pelang eon et al.l (I2008I) . who find a negative corre- 
lation between rarest and z, we do not find any evidence in the 
GBM data for such a correlation. Also, there is no correlation 
between Tgo,Te$t and z when accounting for £p, re st- A Spearman 
rank test for all bursts gives a correlation coefficient of p — 0.04 
with a chance p robability of P = 0.8 1 . However, the p-value in 
Pelangeon et al. is » 10~ 3 which makes it a weak case for such 
a correlation in the first plac e. Our results confirm the analyses 
with Swift detected GRBs (lGreinerll20TTl) . We do note a lack 
of short GRBs for z > 2. This, however, is very likely a selec- 
tion bias, because firstly, short GRBs at such high redshifts must 
be very luminous to be observed by GBM a nd secondly, short 
GRBs are subluminous in the optical band dKann et al J 12008) 
and therefore it is difficult to obtain a redshift measurement. 



4.4. £p, res t vs redshift 

In order to explain the dete ction rate of GRBs at high-z, 
ISalvaterra & Chincarinil (2007) concl ude that high-z GRBs mus t 
be more common (e.g. iDaigne et al.ll2006t IWang & Daill201ll) 
and/or intrinsically m ore luminous (Salvaterra et al. 2009) than 
bursts at low-z (but s ee Butle r et al .I2010h . As already mentioned 
above, Yoneto ku et al.l j2004j) found a tight correlation between 
the 1-s peak-luminosity (L p ) and E p ^ st in GRBs. Assuming that 
the luminosity function of GRBs indeed evolves with redshift 
and that the Yonetoku relation is valid, we would also expect a 
positive correlation of E p resl with z. 
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Fig. 9. Testing the cosmic evolution of 79o, res t. No correlation is 



evident. Bursts with E 



bursts with E 



p,rest 



< 500 keV are shown as crosses, 



/?,rest 



> 500 keV shown as filled circles. The two 



bursts for which are best fit by a simple power-law are shown as 
open diamonds. 



In Fig[10]we present E PtKSl vs z. As was shown in Sect B.l.fl 
and in TableQ] GBM can reliably measure E p down to ~ 15 keV. 
The solid line indicates this redshift-corrected lower limit. The 
Spearman's rank correlation, using only the long GRBs, is p = 
0.58 with a chance probability of P — 2 x 10~ 3 . When including 
the short GRBs, the correlation coefficient effectively remains 
unchanged, whereas the chance probability increases to P — 4 x 
10~ 3 , making a correlation slightly less likely. 

However, this correlation can be explained entirely by se- 
lection effects: GRBs do not populate the empty area in Fig.fTUl 
(low £p jres t and z > 1) because they simply can not be detected 
by GBM. Even though GBM could recover a low E p value of 
such GRBs, as was shown above, the detection of such events 
is very challenging because of the low photon fluxes of these 
events. As one can see in Fig[10]the lower boundary of the ap- 
parent correlation is composed of the bursts that have relatively 
low peak photon fluxes. This is already an indication that these 
events reside at the lower fluence limit for GBM to both trigger 
on these events. 

Since we actually know the intrinsic parameters of the 32 
GRBs, one can test up to which maximum redshift, z ma x these 
bursts could have been detected, i.e. for which GBM would have 
triggered. GBM has many trigger algorithms (various trigger 
time scales for various energy ranges). For the purpose of this 
test, we focus on the 50 keV to 300 range which is the classical 
trigger energy range for a GRB and a timescale of a maximum 
of 4.096 s for long GRBs and 1.024 s for short GRBs. In order 
to shift a GRB to a higher redshift, three observables change: 



1. The duration r 90 (z m ax) = T 90 (z )^ 

) = £ P (zo)i^ 



2. The peak energy E p (z 

3. The flux of the GRB. 



While it is straightforward to account for the changes of 
£p and Tgo, the proper treatment of the flux is more complex. 
RMFIT outputs the spectral parameters, including the normal- 
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10.00 




Fig. 10. Testing the cosmic evolution of £p iIES t for long (dia- 
monds) and short (triangles) GRBs. The solid line indicates the 
redshift corrected lower limit for GBM to measure E p which is 
currently estimated to be ~ 15 keV in the observer frame (see 
Table Q]). Bursts with high (green), intermediate (blue) and low 
(red) peak photon flux (in ph crrT 2 s -1 ) are labeled. 



ization (A^o in ph cirT 2 s keV -1 ) of the spectrum, which can be 
recognized as a proxy for the flux of a GRB. Therefore, in order 
to decrease the flux when shifting the GRB to ever higher red- 
shifts, the normalization has to be decreased accordingly. This 
was done as follows: 

The photon luminosity of a GRB, in the energy band from 
Ei to E2 is defined as follows 



(6) 



L = 47rd 2 L , o N ®i(E)dE, 

j£,/(l+; ) 

where di is the luminosity distance and zo the redshift of 
the burst. The integral describes the photon flux of the burst in 
the Ei to £2 energy range at the rest-frame of the burst. A^o is 
the normalization of the spectrum and <$>i(E) is the shape of 
the spectrum (Band or COMP) with E p - E p (zo), i.e. <5 2 = 
<$>2{E, E p (zo))- The luminosity of a burst is independent of red- 
shift. This, in turn, means that the above equation is valid also 
when exchanging zo with z max . One gets 



■E2/(1+Zmax) 



L = 4nd 



El /(!+--„» 



N Zm ® 2 (E)dE, 



(7) 



where <D 2 = 2 (£, £p(Zmax))- 

Setting the two equations equal and solving for N Zma one is 
left with 



N Zim =N x 



4. „ J* 



E2/(l+zo) 



®i(E,E„(z ))dE 



(8) 



It can be shown, for both the COMP model and the Band 
function, that the fraction of the two integrals is nothing else 
than 



a+ 1 



£l + ^0 2 (E,E p (z max ))dE d+^ +I 



(9) 



where a is the low-energy power law index. 
Thus, 



(i+i 



AT = N x 



4, Zm 



(1 +Zmax)' 
(1 + Z) a+1 



(10) 



For GBM to trigger, two detectors need to be above the trig- 
ger threshold. Therefore, real background information of the sec- 
ond brightest detector of every burst was used. With this back- 
ground data, we shift each bursts in steps of Az = 0.25 to higher 
redshifts and simulate 1000 bursts for each redshift step (chang- 
ing the source lifetime, input E p and normalization as described 
above, additionally adding Poissonian noise to the best-fit pa- 
rameters) by forward folding the photon model through the de- 
tector response matrix, created at the time and location of the 
real GRBs. To determine if GBM would have triggered, we de- 
termine the signal-to-noise ratio (SNR) with 



SNR = 



At • c s - At ■ c b 
VAf • c b 



(11) 



where, Af is the trigger time scale (Af = 4.096 s for long 
GRBs and At = 1.024 s for short GRBs), c s and c h are the 
counts/s of the source and background, respectively. If more than 
90% of the simulated bursts have SNR > 4.5 it can safely be as- 
sumed that the burst would have been detected at this redshift. 
We plot Enrest, which obviously remains constant at all redshifts, 
vs the range of the actually measured to maximum redshift in 
Fig[TT] We already know that the here presented sample of GRBs 
is representative of all bursts detected by GBM (see Sect. 13. 1. Ti l. 
Therefore, we note that the determination of the lowest mea- 
surable E p value is not as crucial as the determination of the 
detector sensitivity. None of the bursts in our sample populates 
the empty region at 2 < z < 8 and 50 < E PtKS t [keV] < 500 
even when shifting them at their maximum detectable redshift. 
We conclude that the correlation between £p, re st and z can be ex- 
plained entirely by the sensitivity limitations of the instrument. 



10.00 




Fig. 11. Same as Fig.[10]but with a line indicating the maximum 
possible redshift, z max that these 30 GRBs can have to still be 
detectable by GBM. 
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4.4.1 . Low number statistics? 

Another way to test if the E p - z correlation is genuine is by 
bearing in mind that it arises simply due to the lack of GRBs 
with measured redshift and E p < 100 keV. In the present red- 
shift sample there are only 4 GRBs with such a low E p value. 
Considering the E p distribution for 488 GBM bursts of the first 
two years of operations (Fig. |2]l it is possible to estimate as to 
how many bursts with E p < 100 keV are expected if 32 GRBs 
are drawn randomly from the full GBM sample. After a run of 
10 4 drawings, a distribution is obtained which peaks at ~ 7. This 
means that, on average, 7 out of 32 GRBs are expected to have 
E p < 100 keV. According to our simulation the probability of 
observing only 4 such GRBs, as we do in our actual sample, 
is P — 8%, which is within the 2<x limit and thus not signifi- 
cant. Therefore, more GRBs with redshift measurement are re- 
quired to understand if the "blank area" at low E p values between 
2 < z < 5 in Fig.[10]is simply underrepresented or not populated 
at all. 



5. Summary & Conclusions 

The Fermi/GBM is a key instrument to study the temporal and 
spectral properties of GRBs. For bursts which have redshift mea- 
surements it becomes even more important since it allows GRBs 
to be studied in their rest-frame. Here, we presented such a study 
for 32 GRBs observed by GBM, focusing on both the temporal 
and spectral properties, as well as on intra-parameter relations 
within these quantities. 

The £p, r e St distribution of the GBM GRB sample covers an 
energy range from tens of keV up to several MeV, peaking at 
as 800 keV. Despite the broader energy coverage of GBM com- 
pared to BATSE, no high-£p population is found. However, the 
GBM E p distribution is strongly biased against XRFs or XRBs 
which have E p values that fall below 15 keV. Additionally, we 
confirm a previously reported parameter reconstruction effect, 
namely that the low-energy power law index a tends to get softer 
when E p is close to the lower end of the detector energy range. 

Using the canonical internal shock model for GRBs, the 
mean £ p , res t (~ 800 keV) implies a dissipation radius of rn as 
1.3 L^ 2 cm. 

Another finding of this work is that the width of the E p dis- 
tribution is, in fact, not c lose to zero. Such a claim has recently 
been brought forward by ICollazzi et al.l (1201 ll) who argue that all 
GRBs are thermostated and thus, have the same £p >res t. If true, 
this would require an unknown physical mechanism that holds 
GRBs at a constant £p ;res t value. However, it could be shown 
here that such a mechanism is not required because our sample 
has an Ep yTest distribution that ranges over several decades in en- 
ergy. 

We find that the rarest stretches from tenths of a second to 
several hundreds of seconds with a median value of m 10 s. 

The Eko distribution ranges from 10 51 erg to 10 55 erg, peak- 
ing at « 10 53 erg. We confirm the E Ptrest - Ei so correlation and find 
a power law index of 0.5, consistent with the values reported in 
the literature but with a significantly larger scatter around the 
best-fit. We also confirm a strong correlation between the 1 s 
peak luminosity of the burst and its E p KSt with a best-fit power 
law index of 0.57. 

We looked for additional correlations between the param- 
eters. We did not find any evidence for a cosmic evolution of 
79o,rest- Although a correlation between Ep jesl and z is not en- 
tirely unexpected from theoretical considerations and looks in- 



triguing on a E p - z plot, we conclude that the apparent relation- 
ship is simply arising due to the detector sensitivity. 
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D. Gruber et al.: Rest-frame properties of 32 gamma-ray bursts observed by the Fermi Gamma-Ray Burst Monitor 
Table 2. Rest-frame parameters of 32 GBM GRBs. 
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Notes. The columns show the name of the GRB, rarest, £p, res t> best fitting model, redshift and redshift references of the 32 GBM GRBs presented 
here. The keyword "Swift" in the mission column denotes whether the GRB triggered both Swift and GBM, whereas "Fermi" indicates that it was 
triggered by GBM only. 

For bursts with * an effective area correction was applied to the BGO with respect to the Nal detectors. 
{a) The r, 0-rest and E pjcst values for GRB 091024 were taken from lGruber et al.l(l201 ll) 

Refe rences. ( 1 ) iTanvir et al l d2010h : (2^ lO'Meara et al.1 d2010h ; (3) ICucchiara & Foxl d20Tol) ; (4) iBergeij dlolob ; (5 ) IWiersema et all d2009bl) 
: (6llCucchiaraet al.U2009aT, : (71 ICucchiara et all (l2009bl) : (8) IXu et all j2009h : (9) leucchiara et al.1 d2009ch : (l O)lLevan et all d2009h : (11) 
Fvnb o et al.N2009[) : (12 ) iMalesani et all (12009177 (13) ICucchiara et al.1 d2009dl) ; (14) ICenko et al l j2009bl) ; (1 5) iThoene et all J20091) ; (16) 
Rossi et all J"2009bh : (1 7) lde Ugarte Postieo et al.1 J2009a|) : (18)lRossi et al.l d2009a|) : (19)lRau et alj 120091) ; (2 0)IChornock et al.l j2009bj) : (21) 
Wiersema et al.ld2009ah ; (22 ) |Tanvir et al.ld2009h ; (23)ICenko et al.l j2009if) ; (24)lchornock et al.1 d2009al) ; (25)lde Ugarte Postigo et al.l( i2009b) 
: (2 6)ICucchiara et al.l d2008h : ( 27)lBerger & Rauchl (l2008h : (28)lD'Avanzo et al J d2008h : (29)lBerger et alj<2008h ; ( 30)ICucchiara & Foxl T2008): 
OD lGreineret al.l (120091) ~(32) lFvnbo et al.l (120081) : (33) [Vreeswiiket all d2008h ; (34) |Prochaska et al.1 d2008T) : (35) lThoene et al.ld2008l) . 
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